function omega = costcompup(sbar,stilde,xlow,xup,p)

r = laborassignment(sbar,stilde,xlow,xup,p);

omega = r.wtop(1)*exp(logcapprod(xup,p))/exp(loglaborprod(stilde,xup,p));

end